Spin pumping by a field-driven domain wall 
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We present the theory of spin pumping by a field-driven domain wall for the situation that spin 
is not fully conserved. We calculate the pumped current in a metallic ferromagnet to first order 
in the time derivative of the magnetization direction. Irrespective of the microscopic details, the 
result can be expressed in terms of the conductivities of the majority and minority electrons and 
the dissipative spin transfer torque parameter /3. The general expression is evaluated for the specific 
case of a field-driven domain wall and for that case depends strongly on the ratio of /3 and the 
Gilbert damping constant. These results may provide an experimental method to determine this 
ratio, which plays a crucial role for current-driven domain-wall motion. 
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I. INTRODUCTION 

Adiabatic quantum pumping of electrons in quantum 
dotsii^ has recently been demonstrated experimentally 
for both charge^ and spir*^. Currently, the activity in this 
field is mostly concentrated on the effects of interactions^, 
dissipation^, and non-adiabaticitj«i. Complementary to 
these developments, the emission of spin current by a 
precessing ferromagnet — called spin pumping — has 
been studied theoretically and experimentally in single- 
domain magnetic nanostructure o^i^d" . One of the differ- 
ences between spin pumping in single-domain ferromag- 
nets and quantum pumping in quantum dots is that in 
the latter the hamiltonian of the electronic quasi-particles 
is manipulated directly, usually by varying the gate volt- 
age of the dot. In the case of ferromagnets, however, 
it is the order parameter — the magnetization direction 
— that is driven by an external (magnetic) field. The 
coupling between the order parameter and the current- 
carrying electrons in turn pumps the spin current^i. The 
opposite effect, i.e., the manipulation of magnetization 
with spin current, is called spin transfer— li^ii^ii^. 

Recently, the possibility of manipulating with cur- 
rent the position of a magnetic domain wall via 
spin transfer torques has attracted a great deal of 
theoreticali^^'i^^'20^'22^^^^^^^^ and 
experimenta]^>2^i^i2ii2^i2&i21i^ interest. Although the 
subject is still controversiali^*^, it is by now established 
that in the long-wavelength limit the equation of motion 
for the magnetization direction $7, which in the absence 
of current describes damped precession around the 
effective field ~SEmm['^]/ (fiSrt), is given by 
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The first is the reactive spin transfer torqu o^^'^^ , which 
corresponds to the term proportional to VfJ on the left- 
hand side of the above equation. It is characterized by 
the velocity Vg that is linear in the curcnt and related to 
the external electric field E by 



Vs = 



\e\Ps 



(2) 



and contains, to lowest order in spatial derivatives of the 
magnetization direction, two contributions due the pres- 
ence of electric current. 



where and denote the conductivities of the major- 
ity and minority electrons, respectively, and Ps is their 
density difference. (The elementary charge is denoted by 
|e|.) The second term in Eq. ([1]) due to the current is the 
dissipative spin transfer torque^ that is proportional to 
^19,20,21 ^ Both this parameter, and the Gilbert damping 
parameter q;g, have their microscopic origin in processes 
in the hamiltonian that break conservation of spin, such 
as spin-orbit interactions. 

It turns out that the phenomenology of current-driven 
domain-wall motion depends crucially on the value of the 
ratio /3/aG- For example, for /3 = the domain wall 
is intrinsically pinnedi^, meaning that there is a criti- 
cal current even in the absence of inhomogeneities. For 
13 /uq = 1 on the other hand, the domain wall moves with 
velocity Vg. Although theoretical studies indicate that 
generically (3 ^ ao^'^i^Si'^, it is not well-understood 
what the relative importance of spin-dependent disorder 
and spin-orbit effects in the bandstructure is, and a pre- 
cise theoretical prediction of P/aQ for a specific mate- 
rial has not been attempted yet. Moreover, the determi- 
nation of the ratio /^/ccq from experiments on current- 
driven domain wall motion has turned out to be hard 
because of extrinsic pinning of the domain and nonzero- 
temperatur o^^i^^ effects. 

In this paper we present the theory of the current 
pumped by a field-driven domain wall for the situation 
that spin is not conserved. In particular, we show that 
a field-driven domain wall in a metallic ferromagnet gen- 
erates a charge current that depends strongly on the ra- 
tio P/aQ. This charge current arises from the fact that 
a time-dependent magnetization generates a spin cur- 
rent, similar to the spin-pumping mechanism proposed 



2 



by Tserkovnyak et al^ for nanostructures containing fer- 
romagnetic elements. Since the symmetry between ma- 
jority and minority electrons is by definition broken in 
a ferromagnet, this spin current necessarily implies a 
charge cmTcnt. In view of this, we prefer to use the term 
"spin pumping" also for the case that spin is not fully 
conserved, and defining the spin current as a conserved 
current is no longer possible. 

The generation of spin and charge currents by a mov- 
ing domain wall via electromotive forces is discussed very 
recently by Barnes and Maekawa^. We note here also 
the work by Ohe et ali^, who consider the case of the 
Rashba model, and the very recent work by Saslow^, 
Yang et aZ4^, and Tserkovnyak and Mecklenburg^. In 
addition to these recent papers, we mention the much 
earlier work by Berger, which discusses the current in- 
duced by a domain wall in terms of an analogue of the 
Josephson effect^. 

Barnes and Maekawa^^ consider the case that spin 
is fully conserved. In this situation it is convenient to 
perform a time and position dependent rotation in spin 
space, such that the spin quantization axis is locally par- 
allel to the magnetization direction. As a result of spin 
conservation, the hamiltonian in this rotated frame con- 
tains now only time-independent scalar and exchange po- 
tential terms. The kinetic-energy term of the hamilto- 
nian, however, will acquire additional contributions that 
have the form of a covariant derivative. Perturbation 
theory in these terms then amounts to performing a gra- 
dient expansion in the magnetization directionii. Hence, 
the fact that Barnes and Maekawa consider the case 
that spin is fully conserved is demonstrated mathemat- 
ically by noting that in Eq. (5) of Ref. [i^ there are 
no time-dependent potential-energy terms. Generaliz- 
ing this approach to the case of spin-dependent disorder 
or spin-orbit coupling turns out to be difficult. Never- 
theless, Kohno and Shibata were able to determine the 
Gilbert damping and dissipative spin transfer torques us- 
ing the above-mentioned method^. Since Barnes and 
MaekawEfc^ consider the situation that spin is fully con- 
served, they effectively are dealing with the case that 
~ P ~ 0- This is because both the Gilbert damping 
parameter ac and the dissipative spin transfer torque pa- 
rameter /3 arise from processes in the microscopic hamil- 
tonian that do not conserve spi n^^i^^i^^'^° . Hence, for the 
case that uq = f3 = our results agree with the results 
of Barnes and Maekawa^. 



The remainder of this paper is organized as follows. In 
Sec. [H] we derive a general expression for the electric cur- 
rent induced by a time-dependent magnetization texture. 
This general expression is then evaluated in Sec. IIIII for 
a simple model of field-driven domain wall motion. We 
end in Sec. IIVI with a short discussion, and present our 
conclusions and outlook. 



II. ELECTRIC CURRENT 

Quite generally, the expectation value of the charge 
current density, defined by j = ~c6H/SA with c the 
speed of light, H the hamiltonian, and A the electromag- 
netic vector potential, is given as a functional derivative 
of the effective action 



(j(x,r)) ^ c 



(SA(x, T 



(3) 



with r the imaginary-time variable that runs from to 
h/{k-QT). (Planck's constant is denoted by h and fceT is 
the thermal energy.) First, we assume that spin is con- 
served meaning that the hamiltonian is invariant under 
rotations in spin space. The part of the effective action 
for the magnetization direction that depends on the elec- 
tromagnetic vector potential is then given byii 

S,S = JdTjd^{jUx,T))A^,{n{K,T))\7o.np{^,T) , 

(4) 

where a summation over Cartesian indices a, a', a" S 
{x, y, z} is implied throughout this paper. In this ex- 
pression. 



j,V(x,r) = — [0t(x,r)r"V„,0(x,r) 
-(V„,</.t(x,T))T"0(x,T)] 

^^^^ A^,ct,H^,r)T''H^,r) , (5) 
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is the spin current, given here in terms of the Grassman 
coherent state spinor = ((/)|, 0p. Furthermore, r" are 
the Pauli matrices, and m is the electron mass. (Note 
that since we are, for the moment, considering the situa- 
tion that spin is conserved there are no problems regard- 
ing the definition of the spin current.) The expectation 
value (• • •) is taken with respect to the current-carrying 
coUinear state of the ferromagnet. Finally, A„(J1) is the 
vector potential of a magnetic monopolc in spin space 
[not to be confused with the electromagnetic vector po- 
tential A(x, t)] that obeys ea,a',a"dAai /dQa" = and 
is well-known from the path-integral formulation for spin 
systems^!. Eq. ([4]) is most easily understood as arising 
from the Berry phase picked up by the spin of the elec- 
trons as they drift adiabatically through a non-coUinear 
magnetization textur o^^'^^ . Variation of this term with 
respect to the magnetization direction gives the reactive 
spin transfer torque in Eq. P]). 

The expectation value of the spin current is given by 

a%(x,r)) = I dr'l dxX,.'(x-x';r-r')^^^^^^ . 

(6) 

The zero-momentum low-frequency part of the response 
function H^^^,(x - x';t - ^') = (j|,"(x, -^ba' (x', r'))^, 
with (• • ■)o the equilibrium expectation value, is deter- 
mined by noting that for the vector potential A(x, r) = 
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— cEe~*"'^/a; the above equation [Eq. ([6])] should in 
the zero-frequency Umit reduce to Ohm's law (js)Q = 
— ft((T| — o-|)E/(2|e|). Using this result together with 
Eqs. (inin]), we find, after a Wick rotation t ^ it to 
real time, that 

(7) 

with V the volume of the system. We note that the time- 
derivative of the Berry phase term is also encountered 
by Barnes and Maekawa in discussing the electromotive 
force in a ferromagnet^. Such Berry phase terms are 
known to occur in adiabatic quantum pumping^. 

We now generalize this result to the situation where 
spin is no longer conserved, for example due to spin- 
orbit interactions or spin-dependent impurity scatter- 
ing. Linearizing around the coUinear state by means of 
n =i {Sn^, Sny, l - Snl/2 ~ 6nl/2) we find that the part 
of the effective action that contains the electromagnetic 
vector potential reads^ 

= JdrJdxJ dr' j dx' J dr" j dx" [Sflaix, r) 

X K„fc(x, x', x"; T, r', r") • A(x", T")<5r!b(x', r')] , (8) 

where a summation over transverse indices a,b £ {x,y} 
is implied. The spin-wave photon interaction vertex 

K,b(x,x',x";r,r',r") = 
A2 



8hc^ 



(9) 



given in terms of the exchange splitting A, is also en- 
countered in a microscopic treatment of spin transfer 
torques^Si. The reactive part of this interaction vertex 
determines the reactive spin transfer torque and, via 
Eqs. ^ and (O, reproduces Eq. ([7]). The zero- frequency 
long-wavelength limit of the dissipative part of the spin- 
wave photon interaction vertex determines the dissipa- 
tive spin transfer torque. (Note that in this approach 
the definition of the spin current does not enter in deter- 
mining the spin transfer torques.) Although Eq. ^ may 
be evaluated for a given microscopic model within some 
approximation scheme^, we need here only that varia- 
tion of the action in Eq. ([8]) reproduces both the reactive 
and dissipative spin torques in Eq. ([1]). The final result 
for the electric current density is then given by 



2\e\V 
d 



(3 / dx^^^^^.V„f2(x,t) 



dt 



dt 



dxAa' (r2(x, t))Va^a'{x, t) 



(10) 



The above equation is essentially the result of a linear- 
response calculation in dfl/dt, and is the central result 
of this paper. We emphasize that the way in which the 



transport coefficients cr| and and the /?-parametcr en- 
ter does not rely on the specific details of the underlying 
microscopic model. Note that the above result reduces 
to that of Barnes and Maekawa (Eq. (9) of Ref. [l^l ) if we 
take (3 = 0. 



III. FIELD-DRIVEN DOMAIN WALL MOTION 

To bring out the qualitative physics, we evaluate 
the result in Eq. pop using a simple model for field- 
driven domain wall motion in a magnetic wire of length 
L. In polar coordinates 9 and 0, defined by n = 
(sin6'cos^, sin0sin0,cos6'), we choose the micromag- 
netic energy functional 

J 



E' 



MM 



K 



Ps 



isin^i 



dx 



' sm 



(vey 



K 



cos + gB cos ( 



(11) 



where J is the spin stiffness, and K± and Kz are 
anisotropy constants larger than zero. The external field 
in the negative z-direction leads to an energy splitting 
2gB > 0. We solve the equation of motion in Eq. ^ 
within the variational ansata^^^ 



(x,t) = 6*0 (x,i) = 2 tan" 



-(rdw(t)-a;)/A 



(12) 



together with 0(x, t) = 4'o{t), that describes a rigid do- 
main wall with width A = \J~TpKl at position rdw(i)- 
The chirality of the domain wall is determined by the 
angle (/)o(i) and the magnetization direction is assumed 
to depend only on x which is taken in the long direction 
of the wire. 

The equations of motion for the variational parameters 
are given b y^^i^^'^^ 

''^dw(i)^ _ 5^ 



A 



A 
K 



aofpoit) = — sin 200 (i) • 



(13) 



Note that the velocity Vg is absent from these equations 
since we consider the generation of electric current by 
a field-driven domain wall. The above equations pro- 
vide a description of the field-driven domain wall and, 
in particular, of Walker breakdown^. That is, for an 
external field smaller than the Walker breakdown field 
= aQKj_/{2g) the domain wall moves with a con- 
stant velocity. For fields B > the domain wall under- 
goes oscillatory motion, which initially makes the average 
velocity smaller. 

Solving the equations of motion results in 



1 



(1 + ' 
gB 



Re 




(14) 
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FIG. 1: Current generated by a field-driven domain wall in 
units of jo — 2L/[\e\{ai — (j|)aG-ft'±], for qg ~ 0.01 and 
various values of /3. The result is plotted as a function of 
magnetic field in units of the Walker breakdown field Bw = 
aGK^/{2g). 



where the indicates taking the time-averaged value. 
Inserting the variational ansatz into Eq. (jlOp leads in first 
instance to 



\e\L 



A 



-Mt) 



(15) 



which, using Eq. (fH|) . becomes 



\ PgB 



1 - 



1 



Re 




2h 



(16) 



As shown in Fig. [1] this result depends strongly on the 
ratio /3/aQ. In particular, for /3 > oq a local maximum 
appears in the current as a function of magnetic field. 
Since ac is determined independently from ferromag- 
netic resonance experiments, measurement of the slope of 
the current for small magnetic fields enables experimental 
determination of /3. We note that within the present ap- 
proximation the current does not depend on the domain 
wall width A. Furthermore, in the limit of zero Gilbert 
damping and /3, the dissipationless limit, we have that the 
current density is equal to (j^) = (cr| — a-^) gB /{\e\L). 
This is the result of Barnes and Maekawa^ that corre- 
sponds to the situation that uq = /3 = 0, as discussed 
in the Introduction. We point out that, within our ap- 
proximation for the description of domain-wall motion, 
putting (3 = ac in Eq. (|16p gives the same result as us- 
ing Eqs. (Uni) and ^ with ac = (3 ^ 0. That the 
situation discussed by Barnes and Maekawa^ is indeed 
that of ac = /? = is seen by comparing their result 
[Eqs. (8) and (9) of Ref. [i^], and the paragraph follow- 
ing Eq. (9)] with our results in Eqs. ((TU]) and (fT^ . 



IV. DISCUSSION AND CONCLUSIONS 



Our result in Eq. (|16p is a simple expression for the 
pumped current as a function of magnetic field for a 
field-driven domain wall. A possible disadvantage in us- 
ing Eq. (|16p . however, is that in deriving this result we 
assumed a specific model to describe the motion of the 
domain wall. This model does in first instance not in- 
clude extrinsic pinning and nonzero temperature. Both 
extrinsic pinningi^ and nonzero temperature^ can be in- 
cluded in the rigid-domain wall description. However, 
it is in some circumstances perhaps more convenient to 
directly use the result in Eq. ([TS]) together with the ex- 
perimental determination of rdw(i)- Since the only way 
in which the parameter (3 enters this equation is as a 
prcfactor of rdw (t) j this should be sufficient to determine 
its value from experiment. We note, however, that the 
precision with which the ratio P/ac can be determined 
depends on how accurately the magnetization dynamics, 
and, in particular, the motion of the domain wall, is im- 
aged experimentally. With respect to this, we note that 
the various curves in Fig. [1] are qualitatively different for 
different values of /3/aQ. In particular, the results for 
f3/a.Q > 1 and (3/aQ < 1 differ substantially, and could 
most likely be experimentally distinguished. In view of 
this discussion, future research will in part be directed 
towards evaluating Eq. (jlOp for more complicated mod- 
els of field-driven domain-wall motion, which will benefit 
the experimental determination of P/aQ. 

A typical current density is estimated as follows. For 
the experiments of Beach et al^ we have that L ~ 20 
^m, and A ~ 20 nm. The domain velocities measured 
in this experiment are fdw ~ 40 — 100 m/s. Taking as 
a typical conductivity (T| ~ 10^ O^^m^^ we find, using 
equation Eq. (|15p with /3 ~ 0.01, typical electric current 
densities of the order of (jx) ^ 10^ — lO"' A m~^. This re- 
sult depends somewhat on the polarization of the electric 
current in the ferromagnetic metal, which we have taken 
equal to 50% — 100% in this rough estimate. Although 
much smaller than typical current densities required to 
move the domain wall via spin transfer torques, electri- 
cal current densities of this order appear to be detectable 
exp er iment ally. 

In conclusion, we have presented a theory of spin 
pumping without spin conservation, and, in particular, 
proposed a way to gain experimental access to the pa- 
rameter (3/aQ that is of great importance for the physics 
of current-driven domain wall motion. We note that the 
mechanism for current generation discussed in this pa- 
per is quite distinct from the generation of eddy cur- 
rents by a moving magnetic domain^. In addition to 
improving upon the model used for describing domain- 
wall motion, we intend to investigate in future work 
whether the damping terms in Eq. ([T]) , or possible higher- 
order terms in frequency and momentum^, have a nat- 
ural interpretation in terms of spin pumping, similar to 
the spin-pumping-enhanced Gilbert damping in single- 
domain ferromagnets^. 
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